In this section, we aimed at exploring the intraspecific variability of tree growth in three hyperdiverse mature forests, encompassing contrasting climates and disturbance regimes across the Tropics. First, for each dataset, we estimated the intraspecific variance of annualized growth using a hierarchical bayesian model. Then, we tested if spatial structure in individual mean growth could be detected. Finally, we tested if intraspecific variability was higher than interspecific variability in individual mean growth.

Note that we here used data on tree growth, which is one component of individual fitness and an imperfect proxy of tree competitive ability, and our analyses could be further tested with other traits in the future.

1 Datasets

1.1 Paracou

The Paracou forest is located in French Guiana. It is one of the best-studied lowland tropical forests in the Guiana Shield region. It belongs to the Caesalpiniaceae facies and has amongst the highest alpha-diversity in the Guiana shield with 150-200 tree species per hectare in inventories of trees with diameter at breast height (DBH) \(\ge\) 10 cm. The studied 93.75-ha area harbours 768 species, 254 genera and 65 families of plants. The Guiana Shield is characterized by Pre-Cambrian granitic and metamorphic geological formations, highly eroded. It is associated with gently undulating landscapes and a very dense hydrographic system. Paracou forest lies in a hilly area, on a formation called “série Armina” characterized by schists and sandstones and locally crossed by veins of pegmatite, aplite and quartz. The topography of the site consists of small hills separated by narrow (< 5 m wide) sandy waterbeds. The altitude varies from 5 to about 45 m above sea level Hérault and Piponiot (2018). The mean annual temperature is 26 °C and winds are generally weak. There is a well marked dry season (from mid-August to mid-November) and a long rain season with a short drier period between March and April (mean annual rainfall of 3,041 mm). Different study programs have been led at the Paracou site (https://paracou.cirad.fr/). Here, we used data from a disturbance experiment, where 15 plots of 6.25 ha were exposed to four different logging intensities between 1986 and 1988. Since then, cartesian coordinates, DBH, species identity and survival of each tree with a DBH \(\ge\) 10 cm have been collected every one or two years, during the dry season (mid-August to mid-November). The final dataset contains ca. 46,000 tree measurements.

Figure 1.1 shows the location of the trees in the Paracou site.

Plots at the Paracou site. Each point is a tree.

Figure S1.1: Plots at the Paracou site. Each point is a tree.

1.2 Uppangala

The Uppangala Permanent Sample Plot (UPSP) is located in South-East Asia, in the Western Ghats of India, and was established in 1989 by the French Institute of Pondicherry in the Kadamakal Reserve Forest, in the Pushpagiri Wildlife Sanctuary, in Karnataka state, India (Pélissier et al. 2011). It is a low altitude (500-600 m) wet evergreen monsoon Dipterocarp forest (Bec et al. 2015). The studied area, of 5.07 ha, is quite steep, with a mean slope angle of about 30–35°. The plots consist in five North–South oriented transects that are 20-m wide, 180- to 370-m long, and 100-m apart center to center, in addition to three rectangular plots which overlap the transects. The transects gather data from 1990 to 2011 and the rectangular plots from 1993 to 2011. The trees with GBH (Girth at Breast Height) \(\ge\) 30 cm (equivalent to ca. 9.5 cm DBH) were measured every 3 to 5 years. The final dataset contains measurements of 3,870 trees of 102 species (including 2 morphospecies). This forest is considered as one of the rare undisturbed tropical forests in the world (Pascal and Pelissier 1996).

Figure 1.2 shows the location of the trees in the Uppangala site.

Plots at Uppangala site. Each point is a tree.

Figure S1.2: Plots at Uppangala site. Each point is a tree.

1.3 BCI

The Barro Colorado Island site is located in central America, in Panama, covered by a lowland tropical moist forest. The zone became an island after a valley was flooded in order to build the Panama Canal, in 1913. It is nowadays considered as the most intensively studied tropical forest in the world. The studied site is a 50 ha plot (500 \(\times\) 1000 m). It has an elevation of 120 m and is quite flat (most slopes are gentler than 10°). Complete censuses of all trees with DBH \(\ge\) 1 cm have been performed every 5 years since 1980. The dataset contains measurements of 328 tree species and 423617 trees, and when only taking into account trees with DBH \(\ge\) 10 cm, it contains measurements of 255 species and 37224 trees. The dataset is available at (Condit et al. 2019).

Figure 1.3 shows the location of the selected trees in the BCI site.

The 50 ha plot of the Barro Colorado Island site. Each point is a tree.

Figure S1.3: The 50 ha plot of the Barro Colorado Island site. Each point is a tree.

2 Data preparation

For the Paracou site, we used the plots 1 to 15 of the Paracou site. We removed all trees of undetermined species and all measures of the years before the perturbations were performed and the biodiversity plots were added (1985-1991).

For the Uppangala site, we removed the census dates which were not common for all plots.

For the BCI site, we removed the trees with missing DBH measures and trees with DBH < 10 cm for consistency with the Paracou and Uppangala sites.

For all three datasets, we computed annualized growth between two censuses as the difference of DBH (Diameter at Breast Height) between two consecutive censuses, divided by the time period between those two censuses. We then removed estimates with a growth < -2 mm/year or > 100 mm/year.

For the Paracou site, after removing 107670 growth estimates for 12629 trees of indeterminate species, this leads to a dataset with 945571 growth estimates for 65805 trees of 614 species.

For the Uppangala site, this leads to a dataset with 57921 growth estimates for 3725 trees of 102 species.

For the BCI site, this leads to a dataset with 167618 growth estimates for 30386 trees of 244 species.

For all three datasets, we finally computed the mean annual growth for each individual tree as the difference of DBH between the first and the last time a tree was measured, divided by the period between those two measurements.

Abundance diagrams for the three study sites. For the Paracou site, only individuals that were identified to the species level are accounted for. The three communities include few dominant species and many rare species (long asymptote).Abundance diagrams for the three study sites. For the Paracou site, only individuals that were identified to the species level are accounted for. The three communities include few dominant species and many rare species (long asymptote).Abundance diagrams for the three study sites. For the Paracou site, only individuals that were identified to the species level are accounted for. The three communities include few dominant species and many rare species (long asymptote).

Figure S2.1: Abundance diagrams for the three study sites. For the Paracou site, only individuals that were identified to the species level are accounted for. The three communities include few dominant species and many rare species (long asymptote).

3 Estimation of IV

To quantify the relative importance of intraspecific and interspecific variability for each site, we used Stan to run a Bayesian hierarchical model, with a species random effect and an individual random effect on the intercept. In order to compare to classic growth models, and estimate an intraspecific variability within diameter classes, we also added an intercept and a diameter fixed effect. We used the brms package (Bürkner 2017, 2018) to implement the model. We scaled the data to improve convergence time.

\(ln(G_{ijt} + 2) = [\beta_0 + b_{0j} + d_{0i}] + \beta_1 \times ln(D_{ijt}) + \epsilon_{ijt},\)

Priors

\(\beta_0 \sim \mathcal{N}(mean=0, var = 1), iid\)

\(\beta_1 \sim \mathcal{N}(mean=0, var = 1), iid\)

\(b_{0j} \sim \mathcal{N}(mean=0, var=V_b), iid\)

\(d_{0j} \sim \mathcal{N}(mean=0, var=V_d), iid\)

\(\epsilon_{ijt} \sim \mathcal{N}(mean=0, var=V), iid\)

Hyperpriors

\(V_b \sim \mathcal{IG}(shape = 10^{-3}, rate = 10^{-3}), iid\)

\(V_d \sim \mathcal{IG}(shape = 10^{-3}, rate = 10^{-3}), iid\)

\(V \sim \mathcal{IG}(shape = 10^{-3}, rate = 10^{-3}), iid\)

Where \(i\) is the individual and \(j\) the species. The model included two fixed effects (\(\beta_0\) and \(\beta_1\)), a random species effect \(b_{0j}\) on the intercept and a random individual effect \(d_{0i}\) on the intercept. \(V_d\) is the intraspecific variance and \(V_b\) is the interspecific variance.

3.1 Paracou

Trace of the posterior estimates for the Paracou site

Figure S3.1: Trace of the posterior estimates for the Paracou site

Density of the posterior estimates for the Paracou site

Figure S3.2: Density of the posterior estimates for the Paracou site

3.2 Uppangala

Trace of the posterior estimates for the Uppangala site

Figure S3.3: Trace of the posterior estimates for the Uppangala site

Density of the posterior estimates for the Uppangala site

Figure S3.4: Density of the posterior estimates for the Uppangala site

3.3 BCI

Trace of the posterior estimates for the BCI site

Figure S3.5: Trace of the posterior estimates for the BCI site

Density of the posterior estimates for the BCI site

Figure S3.6: Density of the posterior estimates for the BCI site

Table S3.1: Variance partitioning of individual growth in three tropical forest sites. Shown are the mean posterior estimates and their estimated error from the growth model fitted for each site
Intercept (\(\beta_0\)) Diameter (\(\beta_1\)) Species variance (\(V_b\)) Individual variance (\(V_d\)) Residual variance (\(V\))
Paracou
Estimate 6.4e-02 3.3e-02 4.7e-01 4.6e-01 7.7e-01
Estimation error 2.4e-02 4e-03 1.7e-02 4.1e-03 2.3e-03
Uppangala
Estimate 9.5e-02 1.9e-01 3.8e-01 6.6e-01 5.9e-01
Estimation error 4.7e-02 1.2e-02 4.4e-02 8.6e-03 1.9e-03
BCI
Estimate 1.8e-01 -2.2e-02 6.7e-01 4.1e-01 8.1e-01
Estimation error 4.3e-02 4.6e-03 3.7e-02 4.1e-03 2e-03
Variance partitioning of individual growth in three tropical forest sites (see Table above).Variance partitioning of individual growth in three tropical forest sites (see Table above).Variance partitioning of individual growth in three tropical forest sites (see Table above).

Figure S3.7: Variance partitioning of individual growth in three tropical forest sites (see Table above).

Overall, a large part of variability in tree growth was imputable to individual effects in the three sites, even after accounting for the effect of diameter on tree growth, showing a high intraspecific variability in growth in this tropical forests (Figure 3.7).

4 Moran’s I analysis

In order to test if individual tree growth was spatially autocorrelated within each site, we performed a Moran’s I test on individual mean growth for each species with more than 5 individuals. For very abundant species (with more than 3000 individuals), we sampled 3000 individuals with a uniform probability. Moran’s I test function was re-written from the ape package’s Moran.I function with a one-tailed test (Paradis and Schliep 2019) in order to be able to select the pairs of neighbours that are less than 100-m apart and thus assess local spatial autocorrelation.

We then computed the proportion of species with a significantly and positively spatial autocorrelated growth (with a 0.05 alpha risk) and the proportion of species with non significant spatial autocorrelation, as well as the corresponding proportion of individuals. The proportions are computed relative to the species on which the test was performed.

Table S4.1: Spatial autocorrelation of the growth of conspecific individuals in three tropical forest sites. Shown are the proportion of species, and of corresponding individuals, in percent, for which individual growth is significantly positively spatially autocorrelated. The spatial autocorrelation of individual growth was tested using Moran’s I index.
Significant Not significant
Paracou
Proportion of species 28.7 71.3
Proportion of individuals 77.7 22.3
Uppangala
Proportion of species 18.5 81.5
Proportion of individuals 45.3 54.7
BCI
Proportion of species 16.9 83.1
Proportion of individuals 51.4 48.6

29%, 18%, and 17% of the species showed a positive and significant spatial autocorrelation in tree growth, in Paracou, Uppangala, and BCI respectively. (Table 4.1).

Although the three forest sites are located in very different areas across the tropics, with contrasting biogeographical history, climate, regime of perturbations and topography, they all showed a substantial spatial autocorrelation of tree growth intraspecific variability, which supports the generality of our results.The magnitude of spatial autocorrelation varied across the three sites however, with Paracou presenting the biggest proportion and Uppangala the lowest, which can be related to each site specificities. The plots of the Paracou site were subjected to disturbance treatments that created numerous large gaps (Molino and Sabatier 2001). This resulted in a high heterogeneity and spatial structuration of the environment. When restricting our analysis to undisturbed control plots (see 6.1 below), the proportion of detected spatial autocorrelation in conspecific growth drops from 78% to 51% individuals, supporting the hypothesis that the more important spatial autocorrelation in Paracou is due to the stronger spatial structuration of the micro-environment associated with disturbance treatments. Conversely, the disturbance regime in BCI is known to produce frequent small gaps (Hubbell 1999), while the mature forest in Uppangala has been reported as barely disturbed (Pélissier et al. 2011). We acknowledge that the hilly topography of Uppangala could make our estimates of distances among individuals based on pairs of tree coordinates provided by field inventory less accurate than in the other two flatter sites.

5 Is variance within conspecifics inferior to variance within heterospecifics ?

In order to test if the performance of conspecifics was locally more similar than that of heterospecifics in the three tropical forests, for each species with more than five individuals and more than five heterospecific neighbours within 100-m radii around these individuals , we computed the semivariances of growth of individuals that are less than 100-m apart, considering conspecifics only on the one hand, and heterospecifics on the other hand. In the first case, semivariance was estimated as the mean of the squared difference in individual mean growth for all pairs of conspecific individuals. In the second case, semivariance was estimated as the mean of the squared difference in individual mean growth for all pairs of individuals with an individual of the focal species and one of another species. For each species, we compared the semivariances obtained with conspecifics only to the ones obtained with heterospecifics using a Mann-Whitney test.

For each site, the proportion of species with a significantly higher or lower semivariance than the one estimated with heterospecifics or for which the difference was not significant, and the corresponding proportion of individuals were then computed.

Table S5.1: Comparison of local intra- and interspecific variability in individual growth for three tropical forest sites. The variability was estimated with the semivariance and the comparison was performed with a Mann-Whitney’s test. The semivariances were computed for all species with > 5 individuals and > 5 heterospecific neighbours within 100 m in the same plot, and considering pairs of individuals that were less than 100 m apart and in the same plot. Shown are the proportion of species, and of corresponding individuals, for which (i) intraspecific variability was significantly lower than interspecific variability, (ii) intraspecific variability was significantly higher than interspecific variability, or (iii) the difference between inter- and intraspecific variabilities was not significant.
Intraspecific variability < interspecific variability (i) Intraspecific variability ~ interspecific variability (ii) Intraspecific variability > interspecific variability (iii)
Paracou
Proportion of species 61.7 39.7 0.333
Proportion of individuals 89.1 10.9 0.0699
Uppangala
Proportion of species 42.2 62.2 4.44
Proportion of individuals 57.7 23.6 18.8
BCI
Proportion of species 46.6 44.3 9.09
Proportion of individuals 75.2 17.2 7.6
Table S5.2: Comparison of local intra- and interspecific variability in individual growth for three tropical forest sites, when controlling for species abundance. See Table 3. Semivariances with heterospecifics were here computed by sampling a maximum of 10 individuals per species.
Intraspecific variability < interspecific variability (i) Intraspecific variability ~ interspecific variability (ii) Intraspecific variability > interspecific variability (iii)
Paracou
Proportion of species 61.7 37 1.33
Proportion of individuals 86.3 12.1 1.56
Uppangala
Proportion of species 42.2 55.6 2.22
Proportion of individuals 61.2 20.4 18.4
BCI
Proportion of species 50 47.7 2.27
Proportion of individuals 71.2 27.8 1.02

In all three tropical forest sites however, growth semivariance was locally significantly higher among heterospecifics than among conspecifics in a large proportion of the tested cases (62%, 42% and 47% of the species, and 89%, 58% and 75% of the individuals in Paracou, Uppangala and BCI respectively; Table 5.1).

Note that, to control for a potential effect of species abundance on the values of semivariances with heterospecifics, we replicate the analysis by computing the semivariances with heterospecifics by sampling a maximum of ten individuals per species. The results were qualitatively unchanged (Table 5.2).

6 Other tests and verifications

6.1 Undisturbed plots in the Paracou site

One of our main hypotheses was that since many environmental variables are spatially structured and that tree growth is largely influenced by environmental variables, tree growth should be spatially structured. Our analysis using Moran’s I test showed that tree growth was indeed spatially structured at our study sites. The Paracou dataset offers the opportunity to test this hypothesis further. Indeed, some plots were disturbed in the early eighties, creating artificial gaps, whereas others were not disturbed. As the creation of gaps results in a strong spatial structure of the light available under the canopy, we hypothesized that growth should be less structured in plots that were not disturbed.

Table S6.1: Spatial autocorrelation of the growth of conspecific individuals in three tropical forest sites. Shown are the proportion of species, and of corresponding individuals, in %, for which individual growth is (i) significantly spatially autocorrelated, (ii) not significantly spatially autocorrelated, or (iii) that have too few individuals to allow a test to be made. The spatial autocorrelation of individual growth was tested using Moran’s I index.
Significant Not significant
Proportion of species 18 82
Proportion of individuals 54.7 45.3

Considering the undisturbed plots only, the proportion of individuals of species with a significant Moran’s I test was 55%, much lower than when including the disturbed plots (89%, see Table 4.1; Table 6.1). This corroborates our hypothesis that the openness of the canopy triggered important spatial structure in disturbed plots due to light gaps, and as in these gaps, trees tend to grow faster, the spatial structure of growth is stronger.

6.2 Smaller stems in the BCI site

In order to be able to compare all three datasets together, we removed all stems that had a DBH inferior to 10 cm in the BCI inventories, which include all stems with DBH \(\ge\) 1 cm. We replicated our analysis on the spatial autocorrelation of tree growth in the complete dataset and found that the spatial structure in individual growth was even more significant.

7 Code implementation

The whole analysis was conducted using the R language (R Core Team 2021) in the Rstudio environment (RStudio Team 2021). The tables were made with the kableExtra package (Zhu 2021), the figures with the package ggplot2 (Wickham 2009), and the code uses other packages of the Tidyverse (Wickham et al. 2019) (tidyr (Wickham 2021), dplyr (Wickham et al. 2021), readr (Wickham and Hester 2020), lubridate (Grolemund and Wickham 2011), glue (Hester 2020), magrittr (Bache and Wickham 2020)) and other R packages (here (Müller 2020), sp Bivand et al. (2013)). Some functions were implemented in C++ thanks to the R packages Rcpp (Eddelbuettel and François 2011, Eddelbuettel 2013, Eddelbuettel and Balamuta 2018) and RcppArmadillo (Eddelbuettel and Sanderson 2014). The pdf and html documents were produced thanks to the R packages rmarkdown Xie et al. (2020), knitr Xie (2014) and bookdown (Xie 2017).

References

Allaire, J., Y. Xie, J. McPherson, J. Luraschi, K. Ushey, A. Atkins, H. Wickham, J. Cheng, W. Chang, and R. Iannone. 2020. Rmarkdown: Dynamic documents for R.
Bache, S. M., and H. Wickham. 2020. Magrittr: A forward-pipe operator for r.
Bec, J. L., B. Courbaud, G. L. Moguédec, and R. Pélissier. 2015. Characterizing tropical tree species growth strategies: Learning from inter-individual variability and scale invariance. PLoS ONE 10.
Bivand, R., E. J. Pebesma, and V. Gómez-Rubio. 2013. Applied spatial data analysis with R. second edition. Springer, New York Heidelberg Dordrecht London.
Bürkner, P.-C. 2017. Brms : An R Package for Bayesian Multilevel Models Using Stan. Journal of Statistical Software 80.
Bürkner, P.-C. 2018. Advanced Bayesian Multilevel Modeling with the R Package brms. The R Journal 10:395.
Condit, R., R. Pérez, S. Aguilar, S. Lao, R. Foster, and S. Hubbell. 2019. Complete data from the Barro Colorado 50-ha plot: 423617 trees, 35 years. Dryad.
Eddelbuettel, D. 2013. Seamless R and C++ integration with Rcpp. Springer, New York.
Eddelbuettel, D., and J. J. Balamuta. 2018. Extending r with C++: A Brief Introduction to Rcpp. The American Statistician 72:28–36.
Eddelbuettel, D., and R. François. 2011. Rcpp : Seamless R and C++ Integration. Journal of Statistical Software 40.
Eddelbuettel, D., and C. Sanderson. 2014. RcppArmadillo: Accelerating R with high-performance C++ linear algebra. Computational Statistics & Data Analysis 71:1054–1063.
Gourlet-Fleury, S., J.-M. Guehl, O. Laroussinie, and E. (Group), editors. 2004. Ecology and management of a neotropical rainforest: Lessons drawn from Paracou, a long-term experimental research site in French Guiana. Elsevier, Paris.
Grolemund, G., and H. Wickham. 2011. Dates and Times Made Easy with lubridate. Journal of Statistical Software 40.
Hester, J. 2020. Glue: Interpreted string literals.
Hérault, B., and C. Piponiot. 2018. Key drivers of ecosystem recovery after disturbance in a neotropical forest: Long-term lessons from the Paracou experiment, French Guiana. Forest Ecosystems 5:2.
Hubbell, S. P. 1999. Light-Gap Disturbances, Recruitment Limitation, and Tree Diversity in a Neotropical Forest. Science 283:554–557.
Molino, J.-F., and D. Sabatier. 2001. Tree Diversity in Tropical Rain Forests: A Validation of the Intermediate Disturbance Hypothesis. Science 294:1702–1704.
Müller, K. 2020. Here: A simpler way to find your files.
Paradis, E., and K. Schliep. 2019. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35:526–528.
Pascal, J.-P., and R. Pelissier. 1996. Structure and Floristic Composition of a Tropical Evergreen Forest in South-West India. Journal of Tropical Ecology 12:191–214.
Pebesma, E. J., and R. S. Bivand. 2005. Classes and methods for spatial data in R. R News 5:9–13.
Pélissier, R., J.-P. Pascal, N. Ayyappan, B. R. Ramesh, S. Aravajy, and S. R. Ramalingam. 2011. Tree demography in an undisturbed Dipterocarp permanent sample plot at Uppangala, Western Ghats of India: Ecological Archives E092-115. Ecology 92:1376–1376.
R Core Team. 2021. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
RStudio Team. 2021. RStudio: Integrated development environment for r. RStudio, PBC, Boston, MA.
Wickham, H. 2009. Ggplot2: Elegant graphics for data analysis. Springer, New York.
Wickham, H. 2021. Tidyr: Tidy messy data.
Wickham, H., M. Averick, J. Bryan, W. Chang, L. McGowan, R. François, G. Grolemund, A. Hayes, L. Henry, J. Hester, M. Kuhn, T. Pedersen, E. Miller, S. Bache, K. Müller, J. Ooms, D. Robinson, D. Seidel, V. Spinu, K. Takahashi, D. Vaughan, C. Wilke, K. Woo, and H. Yutani. 2019. Welcome to the Tidyverse. Journal of Open Source Software 4:1686.
Wickham, H., R. François, L. Henry, and K. Müller. 2021. Dplyr: A grammar of data manipulation.
Wickham, H., and J. Hester. 2020. Readr: Read rectangular text data.
Xie, Y. 2014. Knitr: A Comprehensive Tool for Reproducible Research in R. in V. Stodden, F. Leisch, and R. D. Peng, editors. Implementing reproducible research. CRC Press, Taylor & Francis Group, Boca Raton.
Xie, Y. 2015. Dynamic documents with R and Knitr. Second edition. CRC Press/Taylor & Francis, Boca Raton.
Xie, Y. 2017. Bookdown: Authoring books and technical publications with R Markdown. CRC Press, Boca Raton, FL.
Xie, Y. 2021. Knitr: A general-purpose package for dynamic report generation in r.
Xie, Y., J. J. Allaire, and G. Grolemund. 2019. R Markdown: The definitive guide. CRC Press, Taylor; Francis Group, Boca Raton.
Xie, Y., C. Dervieux, and E. Riederer. 2020. R markdown cookbook. First edition. Taylor; Francis, CRC Press, Boca Raton.
Zhu, H. 2021. kableExtra: Construct complex table with ’kable’ and pipe syntax.